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Abstract. We derive an expression for the accuracy with which sources can be 
localized using a network of gravitational wave detectors. The result is obtained 
via triangulation, using timing accuracies at each detector and is applicable to 
a network with any number of detectors. We use this result to investigate the 
ability of advanced gravitational wave detector networks to accurately localize 
signals from compact binary coalescences. We demonstrate that additional 
detectors can significantly improve localization results and illustrate our findings 
with networks comprised of the advanced Laser Interferometer Gravitational-wave 
Observatory (LIGO), advanced Virgo and Large Cryogenic Gravitational Wave 
Telescope (LCGT). In addition, we evaluate the benefits of relocating one of the 
advanced LIGO detectors to Australia. 



1. Introduction 

The era of advanced gravitational wave detectors [U 121 E] is rapidly approaching, and 
with it the prospect of the regular observation of gravitational wave signals. To fully 
exploit these gravitational wave observations, it will be critically important to correlate 
them whenever possible with electromagnetic counterparts. For example, observation 
of gravitational waves coincident with Gamma Ray Bursts will aid in determining 
the nature of the progenitor [4 j; the gravitational wave signal from core-collapse 
supernovae may carry information about the supernova engine [5]; gravitational and 
electromagnetic observations of binary coalescence will provide independent measures 
of the distance and redshift to the source and consequently enable precision tests 
of cosmology [6 . One method of performing multi- messenger observations will be to 
accurately and rapidly localize the source through gravitational wave observations and 
then use electromagnetic observatories to follow up on the event. Indeed, an ambitious 
project to follow up gravitational wave candidates in a host of electromagnetic 
telescopes is already underway [7 . A single gravitational wave detector gives virtually 
no directional information for a short duration signal. Consequently, localization 
requires a network of detectors and the primary tool for localization is triangulation 
between a signal observed at several sites [8] . 

A number of electromagnetic telescopes are being designed, built and operated 
specifically to search for transient phenomena. These include Palomar Transient 
Factory [9], Pan Starrs [ID] , SkyMapper |TT], LOFAR [12] and many more. These 
observatories will have a field of view of around ten square degrees. There are plausible 
astrophysical mechanisms that produce both gravitational waves and electromagnetic 



Localization with advanced GW network 



2 



signals in various frequencies, for example gamma ray, x-ray, optical and radio (see, 
e.g. [H [33J [14j D2]). Thus, providing a good estimate of the localization capability 
of the gravitational wave network will help inform the follow-up searches which are 
undertaken. In addition, if the source can be localized to a region which contains only 
one, or possibly a few, galaxies then a narrow field telescope could be used to image 
the galaxy in question [T6] . 

In this paper, we investigate the accuracy with which gravitational wave sources 
can be localized using a network of advanced detectors. A network comprising more 
detectors, particularly if they are at widely separated sites, gives the best localization. 
Thus, the recent announcement of funding for the construction of the Japanese LCGT 
[2] is particularly welcome. This detector will augment the advanced LIGO (aLIGO) 
detectors pQ to be installed at sites at Hanford, WA and Livingston, LA in the USA 
and the advanced Virgo [3] detector to be installed in Cascina, Italy and provide a 
four site network. Furthermore, there is significant interest in installing one of the 
three advanced LIGO detectors at the Gingin site in Australia [17J[T8], providing a 
realistic prospect for a five site network of advanced gravitational wave detectors. 

The layout of the paper is as follows. In Section [2] we extend the framework 
presented in [8] to provide a simple, geometrical expression for the localization region 
obtained with any number of detectors. Then, in Section [3J we examine the capacity 
of the various networks to localize sources. The formalism is equally applicable to 
any transient source. However, for simplicity we present a sample of results for binary 
neutron star (BNS) mergers. As noted in [8 , the results for both neutron star-black 
hole and binary black hole mergers with a total mass below around 10 M@ will be 
comparable. 



2. Sky Localization from Triangulation 

For the purposes of this paper, we assume that localization is achieved entirely 
by triangulation of the source by the network of gravitational wave detectors. In 
reality, there is additional information available in the waveform; for instance the 
fact that the gravitational waveform consists of only two polarizations can be used to 
coherently combine the data streams and improve localization p~9j[20]. Indeed, even 
the non-observation of an event in one detector in the network can be used to assist 
with localization. Generically, these effects will serve to improve the localization of 
events. In addition, there will be numerous systematic uncertainties in estimating 
the localization of signals; for example imperfect calibration of the detectors and 
uncertainties in the emitted gravitational waveform. We will not investigate these in 
great detail, but merely note that systematic errors are expected to be smaller than 
statistical uncertainties for all but the loudest sources [8]. 



2.1. Timing 

The accuracy with which the tim^fja gravitational wave signal passed a given detector 
can be determined is given by [21] 



J For an extended signal such as a binary coalescence, this should be regarded as a fiducial time, for 
example the time at which the frequency of the signal passes through 100 Hz. 
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Detector 


Horizon 


Range 


/(Hz) 


a f (Hz) 


a t (ms) at p = 8 


aLIGO NOSRM [23] 


360 


160 


65 


43 


0.46 


aLIGO BNS [23] 


490 


215 


110 


106 


0.19 


Advanced Virgo [3] 


350 


155 


120 


100 


0.20 


LCGT 


365 


160 


100 


88 


0.22 



Table 1. Currently projected advanced detector sensitivity for BNS systems. 
The sensitivity of the detector is encoded in the BNS horizon (distance at which 
an optimally oriented and located signal gives an SNR of 8) and the BNS range 
(the volume and orientation averaged distance at which a BNS gives SNR 8). 
The horizon is a factor of 2.26 larger than the range. The mean frequency / and 
frequency bandwidth give the frequency range in which the detector is sensitive 
to the signal. The frequency bandwidth a f determines the timing width at at a 
given SNR through {T|. The two sets of advanced LIGO numbers correspond to 
no signal recycling (NOSRM) and BNS optimized (BNS) configurations. 

Timing accuracy is inversely proportional to both the signal-to- noise ratio (SNR) p 
and effective bandwidth 07 of the source, defined as 

where the waveform h is normalized such that J \h(f)\ 2 / S(f)df = 1. The 
approximations used to obtain this formula break down at low SNR, where second 
order effects become important [22]. At an SNR of 8, which we consider in several 
later examples, this is about a 10% effect. 

In Table [l] we provide the currently projected sensitivities, frequency bandwidths 
and timing accuracies for the various advanced gravitational wave detectors. 
Numerous configurations have been proposed for the advanced detectors pQ, and it 
is likely that several will be used over the lifetime of the detectors. To illustrate 
the differences, we consider two aLIGO configurations, one with no signal recycling 
mirror (NOSRM) which is likely to be an early configuration, and one optimized 
for BNS detection (BNS) which may be used in later science runs. The anticpated 
noise curves represent the incoherent sum of the principal noise sources as currently 
understood. There will be, in addition, technical noise sources. The sensitivities in 
Table [T] are not the guaranteed performance of advanced detectors, but a good guide 
to the anticipated sensitivity. 

For the most part, advanced detectors will be sensitive to an optimally oriented 
and located BNS to a horizon distance of around 360 Mpc, although the BNS 
optimized aLIGO configuration provides a horizon of almost 500 Mpc. The frequency 
bandwidth of a BNS signal in the detectors will be around 100Hz, leading to a timing 
accuracy of 0.2ms at an SNR of 8. The aLIGO NOSRM configuration has a noise 
curve which rises more sharply at high frequency. Consequently, it has a significantly 
smaller bandwidth and a timing accuracy of almost 0.5ms at SNR 8, more than a 
factor of 2 larger than the other configurations listed. 

For the remainder of this paper, we make the simplifying assumption that all 
the advanced detectors have the same sensitivity and bandwidth. Our standardized 
advanced detector will have a BNS horizon of 360 Mpc (corresponding to a sky and 
orientation averaged BNS range of 160 Mpc at SNR 8) and a bandwidth of 100 Hz. 
It is relatively straightforward to scale the results to other parameter choices using 
Equation 0. 
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2.2. Localization 



The measured time of arrival of a signal in a network of detectors can be used to 
reconstruct the source location. However, consideration of the observed amplitudes 
and phases of the signal at all detectors in the network is critical for extracting the 
full set of parameters, including the distance to the binary and its orientation [24} fT9] . 
For our purposes a gravitational wave signal is described by the location R of the 
source and the time T Q at which the signal passes through the center of the earth. 
Since our primary focus is localization, we take R to be a unit vector describing only 
the position of the source and not its distance. 

The time at which the signal passes through detector i is given by 

Ti = T - R • d, , (3) 

where encodes the separation between detector i and the center of the earth 
(expressed in seconds) . The distribution of the measured arrival times t{ in the various 
detectors, given the actual arrival times T^, is given by 



p(^)=n-Af=T exp 



2n<Ji 



-(U - Ti) 
2a? 



21 



(4) 



where the timing accuracy cr^ for each detector is given by Equation ([!]). 

Measurements of the arrival times in each detector can be used to construct a 
posterior probability distribution for the source's sky location R. This is done by 
applying Bayes' theorem to obtain the posterior distribution for the actual arrival 
times as a function of the observations as 



p(Ti\ti) oc p(Ti)exp 



-(u 



2a 



(5) 



The posterior distribution is the product of the prior distribution p(Ti) for the arrival 
times with the likelihood. Since we are interested in obtaining a distribution for the 
sky location of the event, we would like to re-express ([5| in terms of R. To do so, we 
introduce the measured sky position r and arrival time t Q (in analogy to equation (J3|) 
as: 

ti=t -r'di. (6) 

Strictly, for an event observed in multiple detectors, the set of equations ([6| for t Q and 
r may be over-determined and not allow any solution. However, for a gravitational 
wave signal the observed times in all detectors should be consistent with a common 
geocentric time and sky location. In practice, due to measurement errors, they will 
not be precisely consistent, in which case minimization techniques [25] can be used to 
find the best fit parameters. We will assume that this has been done and proceed to 
eliminate U in favour of r and t Q . Similarly, we make use of ^ to eliminate Ti from 
([5| in favour of R and T Q . The prior distributions are naturally taken to be uniform 
over the sphere (for R) and uniform in time (for T Q ). Finally, after marginalizing over 
T , the posterior distribution for R is 

•J(r-R) T M(r-R) . 



p(R|r) oc p(R) exp 



(7) 



The matrix M, describing the localization accuracy, is given by 

1 



M = 



E 



k a k 



£ 



2a? 



0~A 



(8) 
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where D^- = — dj . A detailed derivation of this result is provided in the Appendix. 
A similar result has been obtained previously in [26] . 

Equation ^ provides a simple extension to an arbitrary number of detectors of 
the two and three detector result given in [8 j. The localization expression has all the 
features we would expect, specifically: localization only depends upon the difference in 
arrival time between the various detectors; localization is improved by extending the 
baseline between detectors and by better timing accuracy in the detectors; the timing 
measurement in a pair of detectors can only serve to restrict the location of the source 
in the direction parallel to the detector separation. The origin of the normalization 
pre-factor for the matrix M arises due to the marginalization over the geocentric 
arrival time T Q based upon timing information at all detectors in the network. 

The matrix M is symmetric and can be diagonalized to obtain three orthogonal 
eigen-directions (e x ,e y ,e z ) with localization accuracies a xl G yi a z respectively. Thus, 
the posterior distribution for the sky location is 



p(R|r) oc p(R) exp 



\x-Xf , (y-Y) 2 




(9) 



where R = (X, Y, Z) are the co-ordinates of the source in our network eigen-basis, 
and r = (x,y,z) describes the measured location. The widths a x , a y , a z encapsulate 
the localization ability of the network. 

The sky position R = (X, Y, Z) is restricted to lie on the unit sphere, and we 
must take this into account to obtain localization regions. Geometrically, equation 
([9| describes ellipsoids of constant likelihood which intersect the unit sphere to give 
the localization distribution as an ellipse on the sky. In most cases, the source will be 
localized to a small enough patch of this sphere that ignoring its curvature is a good 
approximation. In this case, we project M onto directions orthogonal to r using the 
projection 

P(r) =I-rr T , (10) 

where I is the identity matrix. This gives 

M(r) = P(r) M P(r) . (11) 

The eigenvectors and corresponding eigenvalues of M(r) describe the two dimensional 
localization ellipse for the source. The best localization arises when the projection 
preserves two smallest a values, and the worst when it keeps the two largest. In all 
cases, the source is localized with a probability p within an area 

Area(p) * [- ln(l - p)] , (12) 

where o\ and 02 are the localization accuracies of the eigen-directions of the projected 
matrix M(r). 

When the network is degenerate, one or more of the G{ is infinite and the constant 
likelihood ellipsoids become extended. For a three detector network, there is one 
degenerate direction perpendicular to the plane of the detectors and the surfaces 
become cylinders with infinite extent in that direction. Intersecting the cylinder with 
the unit sphere gives two localization regions, one above and the other below the plane 
formed by the detectors. Generically, this degeneracy can be broken by consideration 
of the observed amplitudes in the three detectors, and we assume this can be done. 
Taking the normal to the detector plane to be in the z-direction, the localization area 
for a three site network can be re-expressed in terms of a x and a y as 

Area(p) w 2ira x a y [— ln(l — p)] / cos 9 , (13) 
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where 6 is the angle between the z-direction and the source. If the source lies in (or 
close to) the x — y plane, the it is no longer well localized, and the approximation used 
to obtain ( p~3| ) breaks down. In this case we must go back to the full distribution to 
derive the localization region. For a source in the ^-direction, obtain 

Area(» « 2^a y yfW x [- ln(l - p)\ . (14) 

For a two site network, only a x is finite and the source can only be localized to 
a ring on the sky. The minimal area of the sky containing the source with a fixed 
probability p is independent of the source location and given as: 

V2a x eTr\p) 
Area(p) « 4tt , 

where D is the separation between the two detectors. 

We make repeated use of the localization expressions 
the paper where we examine the localization capabilities of various detector networks. 



12JT5) in the remainder of 



3. Advanced detector networks 

Using the results of Section [2| we can evaluate the sky coverage and localization 
provided by various advanced detector networks. We take the baseline network to 
be the aLIGO network with two detectors in Hanford, WA (denoted H) and one in 
Livingston, LA (denoted L) which we supplement with any number of: 

• the Japanese LCGT detector (denoted J) to be located at the Kamioka mine, 

• the advanced Virgo detector (denoted V) at Cascina, Italy, 

• the installation of an aLIGO detector at the Gingin site in Australia (denoted A), 
and consequently only a single aLIGO detector at Hanford. 

This provides a set of eight advanced detector networks whose sensitivity and 
localization ability can be compared. For simplicity, we take all advanced detectors to 
have a BNS horizon of 360 Mpc and a frequency bandwidth of 100Hz (corresponding 
to a timing accuracy of 0.2ms at SNR 8). We perform three different comparison 
studies. First, in Section |3.1[ we consider a source observed in all detectors in the 
network with an SNR of 8. Then, in Section [3]2j we consider a set of face on binaries 
at a fixed distance from various sky positions. Finally, in Section [331 we consider an 
astrophysically distributed population of sources. A subset of these results were used 
in presenting the scientific case for installing one of the advanced LIGO detectors in 
Australia [171118]. 



3.1. Signal at a fixed SNR 

Table E] summarizes the localization results for BNS sources observed with an SNR of 
8 in all detectors. While it is rather unlikely that a signal will be observed with the 
same SNR in all detectors in a network, this simple scenario is useful for comparing 
different networks. For each of the eight advanced detector networks, we give the 
localization widths (cr x , o y , a z ) from equation ([9| as well as the smallest and largest 
areas of the 90% confidence localization regions. 

The advanced LIGO network (HHL) comprises only two sites and, as expected, it 
is unable to provide satisfactory localization for any sources. The three site networks 
(AHL, HHJL, HHLV) provide good localization in some directions but there remain 
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Network 


Localization^ ) 


Best Area (cleg ) 


Worst Area (cleg ) 


HHL 


(1-4, -) 


1700 


1700 


AHL 


(0.4, 1.7, -) 


8 


150 


HHJL 


(0.5, 2.0, -) 


14 


220 


AHJL 


(0.3, 0.7, 4.8) 


4 


51 


HHLV 


(0.5, 1.4, -) 


10 


150 


AHLV 


(0.4, 0.6, 1.7) 


3 


14 


HHJLV 


(0.5, 0.5, 2.9) 


4 


22 


AHJLV 


(0.3, 0.6, 0.7) 


3 


6 



Table 2. Localization accuracy in networks of detectors. We assume a BNS signal 
observed at SNR 8 in each detector or, equivalently, a timing uncertainty of 0.2 ms 
in each detector. The localization widths correspond to values of (<J X , <J y ,<J z ) from 
equation |9j and have been expressed in degrees. For the two site HHL network, 
localization is only possible in one direction, while for the three site networks there 
is a degenerate direction perpendicular to the plane of the network. T he best and 
worst case areas are 90% confidence areas calculated using equations (|12|15|). 



clear degeneracies which occur for sources close to the plane of the detectors. The 
three site network (AHL) with an Australian detector provides particularly small 
localization regions due to the long baselines between the US and Australia. The 
results for the aLIGO- Virgo (HHLV) network are somewhat better than aLIGO-LCGT 
(HHJL), even though the baselines are comparable; this is because Livingston, Hanford 
and LCGT are not too far from being co-linear, providing one good localization 
direction and one relatively poor one. 

Once a fourth site is added to the network (e.g. AHJL, AHLV, HHJLV), all 
degeneracy is broken. This is most evident in the significant improvement for the 
worst localized sources. However, there remain some directions with relatively poor 
localization for the AHJL and HHJLV networks as these sites are nearly co-planar. 
The five site network AHJLV provides a good localization over the full sky. 



3.2. Source at a fixed distance 

Here, we consider localization for a face on BNS coalescence at a distance of 160 Mpc 
at different locations in the sky. The expected SNR in each detector depends upon 
its sensitivity to the source direction. Consequently, we introduce a simple criterion 
to classify signals as detectable: the expected SNR in the network must exceed 12, 
and the individual SNR at two or more sites must exceed 5. These thresholds are 
motivated by results of LIGO- Virgo searches [27l[28], where an SNR threshold of 5.5 
was used in each detector, and the loudest background events had a combined SNR of 
around 10. For those signals which are classified as detectable, the localization ellipse 
is calculated using the expressions given in Section [2] 

In Figure [I] we plot the localization ellipses for the four and five detector networks 
considered in Section |3.1|§| The results are consistent with those in Section |3T| A three 



site network provides good localization in some regions of the sky, but close to the 
plane of the detectors the localization ellipses become extended. Furthermore, there 
are regions of the sky where the signal is deemed undetectable by our criteria. The 
addition of a fourth site breaks the three site degeneracy and significantly improves 



§ Similar results for a different source model have been presented in |26| . 
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Figure 1. The localization accuracy for face on BNS at 160 MPc in various 
networks of advanced detecors. The ellipses contain the 90% localization regions 
for sources from varioius points in the sky. A x is plotted at points where the 
network would not confidently detect the system. The plots show the localization 
for six different networks: Hanford-Hanford-LCGT-Livingston (HHJL); 
Australia-Hanford-LCGT-Livingston ( AH JL) ; Hanford-Hanford-Livingston- 
Virgo (HHLV); Australia-Hanford-Livingston-Virgo (AHLV); Hanford-Hanford- 
LCGT-Livingston-Virgo (HHJLV) ; Australia-Hanford-LCGT-Livingston-Virgo 
(AHJLV). 

localization, particularly for the worst cases. However, there remain some slightly 
extended regions for the AHJL and HHJLV networks which lie close to the approximate 
plane formed by the four sites. The five site network AHJLV shows excellent resolution 
over the entire sky. 

3.3. A Population of Coalescing Binaries 

The rate of binary coalescences in spiral galaxies is expected to follow star formation 
rate or, equivalently, the blue light luminosity of the galaxies [29], and there may 
well be additional contributions from elliptical galaxies and globular clusters [30] . 
However, with detector sensitivities extending to hundreds of Mpc, it is reasonable 
to assume a population which is uniformly distributed in volume. Additionally, there 
is no reason to expect a preferred orientation of binary systems in the universe, so 
we take a uniform distribution of binary orientations. We use these distributions to 



Localization with advanced GW network 



9 



Network 


Detectable Sources 


Sources Localized within 






ldeg^ 


5 deg 


10 deg 2 


20 deg 2 


HHL 


59 














AHL 


59 


0.4 


5 


13 


30 


HHJL 


85 


0.2 


2 


5 


14 


AHJL 


85 


1 


14 


36 


59 


HHLV 


83 


0.4 


5 


13 


35 


AHLV 


84 


2 


21 


48 


76 


HHJLV 


112 


2 


19 


47 


77 


AHJLV 


114 


3 


34 


84 


111 



Table 3. Sensitivity and localization capability of various different advanced 
detector networks to a population of BNS signals. The number of signals is 
normalized so that there are 40 signals observed in any single detector above 
SNR=8. This corresponds to the annual astrophysical rate estimate presented in 
|29| . For each detected signal, we calculate the localization area and count those 
which are localized within 1, 5, 10 and 20 deg 2 . 



simulate the parameters of a large number of potential signals and for each signal 
determine whether is is "detectable" in a given network, using the same conditions 
as previously. For signals which are detected, we calculate the 90% localization area 
from equations ( T2p5 ). 

Table [3] gives the results of the simulation. For each detector network, the 
expected number of detectable signals as well as the number localizable within 1, 5, 10 
and 20 deg 2 is given. The numbers are normalized to give 40 signals with SNR greater 
than 8 in a single detector, in accordance with the "realistic" estimate of the annual 
astrophysical rate [29]. However, there is at least an order of magnitude uncertainty 
in the rate of BNS signals. Additionally, some relatively simplistic assumptions have 
been made for the detection threshold. Thus, the results in Table [3] should be taken as 
illustrative: significant difference between network performance are meaningful, but 
the actual values should not be taken too literally. 

The results for a population of sources again provide a strong case for constructing 
as many detectors at different sites as possible. As well as an increase in the absolute 
number of observable sources, additional detectors greatly increase the fraction of 
sources which can be well localized. Specifically, for the HHL network no sources 
will be localized within 20 deg 2 . With the introduction of a third site (AHL, HHLV, 
HHJL) a significant fraction (20 to 50%) of sources are localized within 20 deg 2 , and 
the loudest signals may be localized within 5 deg 2 . Moving to a four site network 
(AHLV, AHJL, HHJLV) further improves localization, with the majority of signals 
localized within 20 deg 2 , and as many as 20% to within 5 deg 2 . The five site network 
provides the most remarkable results, with virtually all signals localized within 20 deg 2 , 
a third within 5 deg 2 and the loudest to within a square degree. 

Figure [2] provides a graphical representation of the same data, showing the 
localization distribution areas for different networks. In all cases, the addition of more 
sites improves the localization. Interestingly, for the networks involving an Australian 
detector, the peak of the localization distribution occurs between 5 and 10 deg 2 . This 
corresponds to the typical area of the wide field electromagnetic transient telescopes 
currently being operated or under construction [9] [I2j [lOj [11] . 
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Figure 2. Distributions of localization accuracy for various networks 



4. Discussion 

We have obtained an expression for the localization accuracy of a gravitational wave 
signal in a network of detectors. The localization expression makes use of only timing 
information in the various detectors, with the timing uncertainty taken to be inversely 
proportional to both SNR and signal bandwidth. This extends the results of [8] to 
networks with any number of detectors. As expected, localization only depends upon 
the difference in arrival time between the various detectors, is improved by better 
timing accuracy in individual detector and longer baselines between detectors. The 
expressions presented here make numerous simplifying assumptions by, for example, 
treating the data as Gaussian and stationary, using only timing information for 
localization, taking the leading order contributions to the timing and localization 
distributions, neglecting systematic uncertainties in the waveform and instrumental 
calibration. Thus, while providing useful localization estimates, a real implementation 
of source localization would need to address many additional issues. 

We have examined localization of BNS sources with a network of advanced 
gravitational wave detectors. It seems reasonable to assume that the first operational 
advanced detectors will be at the existing sites (LIGO Hanford, LIGO Livingston, 
Virgo) as much of the infrastructure, including the beam tubes, is already in place. 
Furthermore, there is a good chance that early science runs will be undertaken without 
a signal recycling mirror (the NOSRM configuration in Table]!]). This gives a frequency 
bandwidth of 40 Hz for a BNS signal, compared to around 100 Hz in the other 
configurations. Consequently, the localization accuracy of the network could initially 
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be about a factor of five worse than the results presented in Section [3j meaning only 
the loudest signals will be well localized. However, there is the prospect for significant 
improvement with the installation of signal recycling to provide a more broadband 
spectrum and the addition of new sites (LIGO Australia, LCGT) to the network. 
A four site network operating with broadband sensitivity gives good sensitivity and 
localization accuracy over the majority of the sky, with the potential for as many as 
50% of detected BNS coalescences to be localized (with 90% confidence) within 10 
deg 2 . The five site network provides outstanding coverage over the whole sky with the 
prospect of virtually every signal being localized within 20 deg 2 and a third within 
5 deg 2 . These localization areas are commensurate with the field of view of existing 
and planned electromagnetic transient observatories providing the realistic prospect 
for multi-messenger astronomy in the advanced detector era. 



Acknowledgements 

We have benefited from extensive discussions with many people within the LIGO 
Scientific and Virgo Collaborations. In particular we thank Duncan Brown, Thomas 
Dent, Bangalore Sathyaprakash, Bernard Schutz, Patrick Sutton, John Veitch, Rai 
Weiss and the LIGO Australia committee. We thank Patrick Sutton for independently 
verifying the numbers used in Table [I] This research was made possible thanks to 
support from the Royal Society. 



Appendix A. Detailed derivation of the localization expression 



Here, we provide a detailed derivation of the localization expression Q presented in 
Section 2.2 We begin with the probability distribution for the signal arrival times Ti 
at the different detectors, given the measured times ti ([5|, 



p(Ti\ti) cxp(T,)exp 



2a? 



(A.l) 



In order to proceed, we assume that the measured times U are consistent with a signal 
from a sky location r and geocentric arrival time t Q . For an event observed in multiple 



detectors, the set of equations (A.2) for t Q and r may be over-determined and not 



allow any solution. However, for gravitational wave signals the observed times in all 



detectors should be consistent. In this case, we can replace the ti terms in ( A.l ) using 



to-r-di. 



(A.2) 



Similarly, we can replace the actual arrival times Ti with the actual sky location R 
and geocentric arrival time T using 



R-di 



(A.3) 



Since the source must have originated from a specific point on the sky, there is no 
complication in doing so. 

Having eliminated both the measured and actual arrival times in favour of the 
sky location and geocentric time, the posterior distribution (A.l) can be re-expressed 
as 



p(H,T \r,t o ) (xp(R,T )exp 



[(t - T ) - (r - R) ■ d,] 2 
2a? 



(A.4) 
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Next, we would like to marginalize over the arrival time T^. The measurement 
uncertainties in the arrival time at the various detectors are typically fractions of a 
millisecond. Even in cases where the time of the event is well known from another 
astronomical observation, the arrival time of the gravitational wave signal will not be 
known with millisecond accuracy. Consequently, it is appropriate to take a uniform 
prior on the arrival time T , and also make the priors on arrival time and sky location 
independent. Then, the distribution (A.4) can be written as 



p(R,T |r,t o ) ocp(R)p(T ) exp 



(to 



2(t -T )J2 



(r - R) • di 



(A.5) 



[(r-R)-^ 



2af ^ 2af 

In order to marginalize over T , we simply integrate over that variable. This is most 
easily done by completing the square in the expression (A.5), and performing the 

(A.6) 

r / -1 \ 

oc p(R) exp 



Gaussian integral. The resulting distribution for the sky location R is 
p(R|r,t ) = J dT o p(R,T \r,t o ) 

-1 



[( r -R).d,] 2 ^ 



■ p(R) exp 



--(r-R) T M(r-R) 



£ 



(r — R) • di 

2a? 



-I 2N 



where we have introduced 



M : 



1 



r -2Zj 9/t 2 /t 2 and D 



(A.7) 



This is the form which is given in the Equations ^ and ([8| in Section 2.2 and used 
throught the paper. 
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